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Abstract 

We derive the kinetic equation for fermions and antifermions interacting with a planar Higgs 
bubble wall during the electroweak phase transition using the 'evenisation' procedure. Equations 
of motion in a relativistic quantum theory do not mirror classical relations unless one uses evenised 
operators. We give a brief introduction to evenisation and then use the evenised Heisenberg 
equations of motion to obtain the velocity and force for the particles in the presence of the Higgs 
bubble wall. Keeping quantum contributions to 0{h) in the equations of motion we obtain the 
semi-classical force obtained earlier by other techniques. 
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I. INTRODUCTION 



Electroweak baryogenesis allows the possibility that the observed matter-antimatter 
asymmetry of the universe may have been obtained during the electroweak phase transition 
at temperatures of around 100 GeV, corresponding to an energy scale that is experimentally 
accessible. Hence there has been considerable interest in ascertaining the details of elec- 
troweak baryogenesis models, in particular, the kinetic equation or Boltzmann equation for 
the transport of particles through the Higgs bubble wall in a first order electroweak phase 
transition. 

It is well known that the characteristics of the Boltzmann equation, in the absence of 
collisions, are the single particle trajectory equations, i.e. if /(p, x, t) is a distribution 
function for a system of particles then '^'^^^f^'*'' = 0, where x and p satisfy single particle 
equations of motion. Therefore expanding the total derivative using the chain rule as 

dtf + ^d^f + ^dpf = 0, (1) 

the quantities ^ and ^ can be obtained for the system from the single particle equations of 
motion (Hamilton's equations). However, in relativistic quantum theories it is non-trivial to 
obtain the single particle equations of motion from the Hamiltonian because of interference 
between particle and antiparticle states causing Zitterbewegung. Obtaining equations of 
motion that can have a consistent single particle interpretation requires the use of 'evenised' 
operators |l| , as we discuss below. We apply the procedure of evenisation to calculate particle 
trajectories to 0{h) in the presence of the Higgs bubble wall during the electroweak phase 
transition (presumed to be first order). We then use this to obtain the kinetic equation. 



In Ref. P] the energy relation, u;s±(p,x), and a kinetic equation in the absence of colli- 
sions for particles and antiparticles interacting with the Higgs bubble wall are obtained from 
the equation of motion for the Wigner function G^{x, A; ), which is the Wigner transform of 
the Wightman function i (^p{x')'ijj{x")). (Also see Ref. 3|.) In this approach, particle and 



antiparticle distribution functions are extracted from the Wigner function using a spectral 
decomposition that separates positive and negative energy states. The particle and antipar- 
ticle distribution functions are then found to satisfy a coUisionless quantum corrected kinetic 
equation in the presence of the external force field provided by the Higgs wall. 

The kinetic equation has also been derived in Ref. Q| . The energy relation was obtained 
using an WKB ansatz for components of the (spinor) wave function. The single particle 



trajectory was obtained by treating the energy relation as the Hamiltonian. Earher analyses 
of the kinetic equation using the WKB ansatz were developed in Refs. jsi 1^ . 

In the above approaches the dominant contribution to the CP violating source that 
leads to a baryon asymmetry is associated with the force term in the kinetic equation. 
Alternate formalisms for the derivation of the kinetic equation have been provided in Refs. 
QBSHInl. In Refs. and a source term was included in the diffusion equation to 



describe the CP violating interaction between the fermions/antifermions and the bubble wall 
in the presence of the thermal plasma. In the former, the source was calculated by including 
reflection and transmission of particles and antiparticles by each layer of the wall. In Ref. 
^ the expectation value of the CP violating current was calculated using the Closed Time 
Path formalism and its divergence was included as a source. In Refs. 0, Q| the Closed Time 
Path formalism was used to derive the kinetic equation itself from the equation of motion 
of the relevant Green's function. 

The kinetic equation obtained below by us agrees with Ref. j^. In particular we obtain 
the semi-classical force term, first discussed in Ref. Q]. Our approach provides an alter- 
nate and intuitive derivation of the (collisionless) kinetic equation relevant to electroweak 
baryogenesis to 0{h). Furthermore, we are not aware of earlier attempts to use evenisation 
to study particle trajectories at higher than 0{h^). 



II. 



A first order electroweak phase transition proceeds via the formation of Higgs bubbles. 
As the bubbles expand they move through the ambient sea of quarks, leptons and other 
particles. To calculate the baryon asymmetry created as the ambient particles interact with 
the expanding Higgs bubble it is necessary to study the kinetic equation describing the 
passage of the particles through the bubble wall. Below we first obtain an expression for 
the single particle energy to 0{h) in terms of the evenised position and momentum. We 
then obtain expressions for the evenised velocity and force and substitute these in Eq. ([T} 
to obtain the kinetic equation. 

We first give a brief introduction to evenisation. In relativistic quantum mechanics the 
eigenvalues of many operators do not agree with their classical values. For example, for a 
free particle the eigenvalues of the velocity operator, d'k/dt = {ih)~^[ii., H] = ca, are ±c 
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even for a massive particle. In addition, one can not obtain relations analogous to classical 
equations for expectation values of quantum operators (i.e. Ehrenfest's theorem), as also 
evidenced by the expression for the velocity operator (quantum mechanically d{x.) /dt = c{a) 
while classically, d:K/dt = p/E). All this is because of interference between particle and 
antiparticle states leading to Zitterbewegung One may resolve this by considering 'even' 
parts of operators, and by replacing operators by their 'even' counterparts in equations of 
motion. This then allows a consistent single particle interpretation of the Dirac equation. 

As an illustration, the evenised velocity operator is given by [dii/dt] = [a] = 
pA/i/p^ + m^, where A is a sign operator with eigenvalues ±1 (and we now set c = 1). This 
mimics the classical relation. Analogously, we shall obtain an expression for the Hamiltonian 
(actually its square) in terms of evenised position and momentum operators and use that to 
obtain an expression for the energy of the particle. We will further obtain expressions for 
the evenised velocity and force on the particle in the presence of the bubble wall and use 
these to obtain the kinetic equation. 

Any operator A can be split into an even part A and into an odd part |^}- Even and 
odd parts of the operators are defined by using the sign operator 

A = ^, (2) 



where H is the Hamiltonian. The eigenvalues of the sign operator are ±1, corresponding to 
particle and antiparticle states. Then 

A 



i (i + AiA) (3) 

{A} = UA- aAa) = 1[A,a]a (4) 



2 V / 2 

The even part of the product of two operators A and B can be written as 



AB 



A 



B 



+ {a} {b} (5) 



We now apply the above to the system under consideration. 

The Lagrangian describing the interaction of particles with the bubble wall can be mod- 
eled by 

C = iip pip — mRi/ji/j — imiip'j^ip (6) 
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where mn and m/ are real numbers.^ The Higgs bubble will be treated as a background field 
which provides a complex spatially varying mass for the particles in the bubble wall frame, 
in the limit of large bubbles when the walls can be treated as planar. The corresponding 
Dirac equation 

(i p -niR - imi'-f^)^) = (7) 



can be rewritten in the form 



.^^H,. (S, 



with H given by 

H = a ■ p + (^TTiji + i7°7^m7 . (9) 

We now use the method of evenisation and evenise H^. We discuss later why we evenise 
rather than H. 

We first write after some algebra, 

^2 ^ p2 + |^|2 +a(3[-p,mji] + iaf^^[-p,mj] (10) 

where, |mp = m|. + m|. Note that a and p are contracted in the last two terms on the 
right hand side. We take the wall to be planar in the x — y plane and so m is a function of 
z only. Then the Hamiltonian can be written as 

= p"^ + Imp - {a^f)ff{-ihd,mji) -i{a^f)f3{-ihd^mj), (11) 

where z = and = d/dx^. Defining the spin operator S as a^^ (without the standard 
h/2 in the definition to be consistent with the notation in Refs. ; the spin operator 

defined here is S with Pauli matrices along the diagonal) the Hamiltonian can be rewritten 
as 

= + |m|2 - S^^^f{-ihd,mR) -iS^f3{-ind,mj) . (12) 

We chose to express a^(3 and d^7°7^ in terms of 5*^ as this simplifies the evenisation of the 
corresponding terms below. 

We now evenise [H"^] to order h. The spin operator in general does not commute with 
the Hamiltonian. In order to simplify the problem, as in Ref.Q, we assume that we are 



^ Alternate ways of writing the mass terms are — jmlV^e'^^^V' ^^^d —m'ipP]^'ip — m*ipPLip, where m = \m\e'' 
ran + imi, as in Refs. |2, 
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in an inertial reference frame where x and y components of the momentum are zero and 
consequently and will be absent from the Hamiltonian. This will allow us to set {S^} 
to later. As mentioned earlier, we wish to express [H"^] in terms of [pz] and [z] (where 
Pz = p^). The odd part of Pz will be proportional to [A,^^] (see Eq. above) and so 
will be of order h. Therefore, using Eq. (0), we can approximate [p^] by [pzY'- Similarly 
we approximate [m\{z) + mj(i)] as [mii{z)Y + [mi{z)Y'. Furthermore, since {z} is 0{h), 
[m/j(i)] = mjidz]) + 0{h^) as any dependence of the even operator [mR{z)] on the odd 
operator {z} can only appear as {z}^- 

The last two terms on the right hand side of Eq. ()12j) are 0{h) and therfore we will use 
the sign operator defined upto 0{fi^) to evenise these terms. We define the zeroth order 
energy as 



Eo = ^Jpz^ + |m|2 (13) 

where and |mp are real numbers and are expectation values of the corresponding oper- 
ators in an eigenstate of definite energy and spin, ("designates an operator.) With this we 
define the sign operator upto O (ji^^ as follows 

^» = / I. (") 

/Pz + 



Note that (Aq)^ = l + 0{h) which we shall use in the derivation of evenised operators below. 
Now 

[S'] = [a^f] = , (15) 

i.e., {S^} = as 5*^ commutes with the Hamiltonian (which is true only in this preferred 
inertial frame), and 

Eo 



[ff] = -^^Ao + 0{h) (16) 



0] = 9^Ao + Oih) . (17) 
Derivations of the above expressions are provided in the Appendix. Thus 

[H'] = [pz? + rnl + m'j- hS'^^dziunk, - hS'^^dzmjAo (18) 

where all mass functions are functions of [z]. Replacing evenised operators [z] and [pz] by 
real numbers representing their expectation values in states of definite energy and spin. 
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\E, s) for particles and | — E, —s) for antiparticles, we deduce the corresponding expression 
for the energy to be 

hs 

E'^ = + |m|^ - —{mjd^mR - niRd^mi) , (19) 

-C/O 

where we have replaced by its eigenvalues ±s and Aq by ±1 for particle/antiparticle 
states. Writing itir = \m\ cos6 and m/ = |m| sin^ one can rewrite the above equation as 

= p/ + |m|2-^|m|2^' (20) 

Eq 

The energy relation for both particles and antiparticles can then be written upto order h as 

E = ^p,2+|^|2_^|^|2^/ (21) 

This expression for the energy in terms of the position, momentum and spin of the parti- 
cles/antiparticles is the same as the energy relation obtained in Ref. P]. Note that if we 
had wished to obtain the energy relation to 0{h) by evenising H instead of evenizing we 
would have faced a problem as we would have needed A to 0{h) which itself requires E to 
0{h). However, since {H} = 0, [H] = [H"^]^^"^. We are now able to define the sign operator 
upto 0{h) as A = H/E and use the same to obtain the evenised velocity and force to 0{h) 
for the kinetic equation. 

We now obtain the kinetic equation. Using the chain rule for partial differentiation the 
kinetic equation, in the absence of collisions, is written as 

dtfs± + ^djs^ + ^d,Js± = , (22) 

where fs± are the particle and antiparticle distribution functions. We will associate dz/dt 
and dpz/dt with the expectation values of [dz/dt] and [dpz/dt] in states of definite energy 
and spin as before. Implicitly we are assuming here that the form of the quantum Boltzmann 
equation is the same as that of the classical Boltzmann equation and the quantum corrections 
are contained only in the coefficients of the equation, namely, in the expression for the 
velocity and the force. The evenised expressions for the velocity and force to 0{h) are 
derived in the Appendix as 

[dz/dt] = [~^[^,m] 



\mf , hS'Aoi\mm:\j^^ (34) 



V 2E 2^2 

where \m\ and 9' are functions of [z] and E is as given in Eq. (PT|) . Substituting the 
eigenvalues of the above operators in the kinetic equation 2, we get 

a,/.. + + (-^ + ^^^^) ^.U = . (25) 

We note that this agrees with the kinetic equation obtained in Ref. P|. ^ 

Several comments are in order here. If we apply Hamilton's equations to the energy rela- 
tion in Eq. (j2H) . i.e., dz/dt = dE/dpz and dp^/dt = ~dE/dz, we get expressions for dz/dt 
and dpz/dt that differ from the expectation values of the corresponding evenised operators 
[dz/dt] and [dpz/dt]. There can be two reasons for this. Firstly, the real number pz appear- 
ing above in Eq. ()21|) is not the canonical momentum. It represents the expectation value 
of [pz] whereas the canonical momentum would be associated with the expectation value of 
Pz- Secondly, ignoring for now the distinction between pz and [pz], Ehrenfest's theorem for 
a relativistic system would imply d{z)/dt = {dH/dpz) and d{pz)/dt = —{dH/dz), and not 
d{z)/dt = ^{H)/^{p^ and d{pz)/dt = -d{H)/d{z). 

In the literature 2, 0] there has been discussion on identifying the canonical and kinetic 
momentum of the particle. One may presume that since pz is identified with the —idz 
operator it is conjugate to z and represents the canonical momentum. Now our energy 
relation agrees with that of Ref. P]. But the energy relation in Ref. S is necessarily in 
terms of the kinetic momentum as the Wigner function invoked in Ref. l3| is always written 
in terms of the kinetic momentum of the system (see Sec. 3 of Ref. |l3|). This indicates 
that our real number pz = {[pz]) in the energy relation above is the kinetic momentum of 
Ref. |3|j|. It is not obvious to us why the evenised canonical momentum in this problem 
becomes equivalent to the kinetic momentum. 



^ The overall minus sign in the expectation value for the velocity and force for antiparticles can be absorbed 
by taking the expectation value of the momentum in antiparticle states to be —pz , as for systems with 
constant mass, where the energy eignenstates are also momentum eigenstates. 

^ In the earlier version of this paper, we had identified antiparticles of spin s with negative energy particle 
solutions of spin s, rather than —s (see Sec. 7.1 of Ref. Hence the expressions for energy and force 
were different for particles and antiparticles of the same spin. Our earlier results agreed with Refs. 0,0 
which perhaps contain the same misinterpretation. Our current results reflect the P and CP violation 
and C conservation properties of the Lagrangian in Eq. © . 
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If one wishes to define a classical Hamiltonian function which gives the expressions for 
velocity and force above by applying Hamilton's equations, one may define a canonical 
momentum by the following ansatz of Eq. (61) of Ref. 14 1 

hse' \ 



Pc 



Pz\l + 



2{pI + |m|2) 

and reexpress the energy relation in terms of pc as 



(26) 



E 



; (1 H ) + \mr 



hs9' 



Ip, 



Pc + l^^l ~ 



hs9' 



(27) 
(28) 



Then dz/dt = dE/dpc and, rewriting Eq. aspz = Facpc, dpz/dt = {dFa.c/dt)*pc — Fac* 
dE/dz. These agree with the expressions of the evenised velocity and force above. Eq. ()26|) 
agrees with the relation between kinetic and canonical momentum in Ref. ^, with acp of 
Ref. 1^] set to 0. 



III. CONCLUSION 



In conclusion, we have provided an alternate procedure for obtaining the energy relation 
and the kinetic equation for fermions and antifermions interacting with the Higgs bubble 
wall during the electroweak phase transition using the method of evenisation. Our derivation 
of the single particle behaviour using the Heisenberg equations of motion has allowed us to 
include quantum corrections systematically to 0{h). This provides a direct and intuitive 
way of obtaining the kinetic equation in this limit. Our results agree with those that may be 
obtained using more formal quantum field theoretic methods involving the Wigner function, 
or the WKB method. In particular, we rederive the semi-classical force that may be obtained 
by these other methods. 
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Appendix 

We present below the derivations of the various evenised expressions quoted in the text. 
Standard formulae that we shall use are = — (7°7^)^ = (a^)^ = 1- The symbols square 
brackets are used below for both commutators and evenised operators but the difference is 
clear from the context. 

Evenised /3 and to 0(^°): 



2[/?] = /? + Ao/3Ao ^ 

^ q; V + Mr + il^l^rhi ^ <- 

— PH ^ P-'i-o 

_ « , « -« V + - ij^j^rhi - 

— P + P ^ -'i-O 

— P + P ^ -'^0 

= ^A., ° (29) 

where we have used Aq = 1 + O(^) in the last equality. Since the Ihs is even and the sign 
operator on the rhs is even, mji{z) may be replaced by [mfl(z)], which, as explained in the 
text, is equal to mij([z]) to 0{h). Therefore 

m = ^^^Ao + 0{h) . (30) 

Similarly 

2[7°7'] = ff + koffko 

= 7 7 +7 7 Ao 

= -^A.. (31) 

Therefore, replacing mi{z) by m/([i]), 

[T»fl = -^A, + 0(R). (32) 

Evenised expressions for the velocity and the force to 0{h): 
Since we wish to work till 0{h) we now use the sign operator A defined as 

A = § (33) 



with E defined to 0(h) as in Eq. Now dz/dt = -{i/h)[z,H] = a^. Then 

E 

~ E 



. .~AE + 2a^p\ 
+ —A 



-fA^fA, (34) 

Since the Ihs of the above equation and the sign operator on the rhs are even we may replace 
Pz by [pz]- Therefore we get 

[dz/dt] = %iA . (35) 
E 

The force dpz/dt = —{i/h)[pz, H] = —(3m'f^ — i^^^^m'j. Then 



2[/3my = (5m'ji + Apm'j^A 



^ ^-A . (36) 

Again, we may replace •m'^iriji by [m'j^rh^] as both sides of the equation should be even. As 
{m^jlm^} is 0{h'^), Eq. (0) implies that this reduces to [m'j^{z)][mji{z)]. This may further 
be rewritten as m'j^{[z])mj^{[z]) , since as discussed earlier in the text, the dependence on 
{z} of each term in the product will be at 0{fi?). For the second term in the numerator 
we replace (3a^m'^ by [fSa'^rh'^. Since this term is already of 0{h) from the commutation 
employed in its derivation, we can then write the product as [m'^(i)] = [/5d^]m'^([i]). 
Now = S'^ff and since {S^} = 0, 0a^] = S^[^^^^] = -zS^^^^Aq. Therefore 

\ 2E 2E J 

For the second term in dpz/dt, 

2[ffm'j] = ffm'j + Affm'jA 

\e + 



YYm'j + Yl rn'j — — A - ^ 



E 

7°7^d^ = Invoking arguments as above and substituting for [[3] we get 

^fr^^ = MM. _ ^S^Ao^MWMM) A . (39) 
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Thus 

[dp./dt] = (^-^ + ' A , (40) 

where we have used rhjirh'j — rhirh'^ = (|mp^')'. |m| and 9' above are functions of [z]. 
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